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Abstract 

In this work we study the Lepton Flavour Violating semileptonic r decays: 1) 
r -» uPP with PP = ir + ir-,ir°ir ,K + K-,K K ; 2) r pP with P = ir°,n,n' and 
3) t — ► fiV with V = p°,4>. We work within the context of two constrained MSSM 
scenarios: the CMSSM-seesaw and NUHM-seesaw, with a MSSM spectrum extended 
by three vr and their SUSY partners and where the seesaw mechanism for neutrino 
mass generation is implemented. A full SUSY one-loop computation is presented and 
the importance of the various contributions, the 7-, Z-, and Higgs bosons mediated 
ones, are analysed. The hadronisation of quark bilinears is performed within the chiral 
framework. Some discrepancies in the predicted rates for BR(r — > ui]), BR(r — > fin') 
and BR(r — * fj,K + K~) are found with respect to previous estimates, which will be 
commented here. These three channels will be shown to be the most competitive ones 
to test simultaneously the Lepton Flavour Violation and the Higgs sector. We further 
present here a set of approximate formulas for all the semileptonic channels which we 
believe can be useful for further comparison with present and future data. 
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1 Introduction 



Lepton Flavour Violating (LFV) processes provide one of the most challenging tests of su- 
persymmetric (SUSY) extensions of the Standard Model (SM) of Particle Physics [1-4]. One 
of the most popular ones among these extensions is the Minimal Supersymmetric Standard 
Model (MSSM) enlarged with three right handed neutrinos and their corresponding SUSY 
partners, and where the physical neutrino masses are generated via a seesaw mechanism [5, 6]. 
Within this SUSY-seesaw context, the light neutrino masses and neutrino mixing angles can 
be easily accommodated in agreement with present data [7] by setting appropriate input val- 
ues for the heavy right handed neutrino masses, within the range Mr ~ (10 10 — 10 15 ) GeV, 
and appropriate Yukawa couplings, Y v . The hypothesis of Majorana neutrinos is crucial in 
this concern, because it is only for them that large Yukawa couplings, say Y v ~ 0(1), can be 
set. An interesting connection between neutrino and LFV physics then follows, because the 
large Yukawa couplings induce, via loops of SUSY particles [8], important contributions to 
the rare LFV processes. In fact, these contributions are in some cases, already at the reach of 
their present experimental sensitivity. So far, the most sensitive LFV process to the Yukawa 
couplings in the SUSY-seesaw context is \x — > cy, where the present experimental sensitivity 
is at 1.2 x 10~ n [9,10]. In the future, if the announced improvement in the sensitivity of 
/j — e conversion in nuclei of up to 10 -18 is finally reached [11], this process will be by far the 
most competitive one. Regarding the tests of Lepton Flavour Violation (LFV) in the r — /i 
sector, the most competitive one at present is r — > /ry, whose upper bound is now set to 
1.6 x 1CT 8 [12-15]. Furthermore, the sensitivity to LFV in r — > 3/i has also improved notably 
in the last years. The present upper bounds from BELLE and BABAR collaborations are 
3.2 x 1CT 8 [16] and 5.3 x 1CT 8 [17], respectively. This leptonic channel has the advantage 
over the radiative r — > 7x7 decay that provides a test not only of SUSY but also of the 
Higgs sector. It is remarkable that both r — > 3/i decay [18-22] and fx — e conversion [23, 24] 
in nuclei can get important contributions from Higgs mediated diagrams in SUSY scenarios 
with large tan (3 and light MSSM Higgs bosons. 

In the present work, we study the LFV semileptonic tau decay channels which are also 
of interest because of the recently reported sensitivity by BELLE and BABAR collabora- 
tions [25-28] that are, for some channels, already competitive with the LFV tau leptonic 
ones. In particular we analyse here the following semileptonic tau decays: 1) r — > fiPP 
with PP = 7T+7T-, 7i°7i°, K + K-,K°K°; 2) r -> fiP with P = TT°,r],r]' and 3) r -> fjV with 
V = p°,0. Their present upper experimental bounds (90% CL) are summarised in Table 1. 
We perform a full one-loop computation of the rates for all these processes within the context 
of two constrained SUSY-seesaw scenarios which are of particular interest: The usual con- 
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strained MSSM-seesaw (CMSSM-seesaw) scenario [29], with universal soft SUSY masses at 
the gauge coupling unification scale, and the so-called Non-Universal Higgs Mass (NUHM) 
scenario [30], with all those soft masses being universal except for the Higgs sector ones. In 
this later case the predicted Higgs particle masses can be low, indeed close to their present 
experimental lower bounds (for the SM Higgs the present bound is m H > 114.4 GeV 95% 
C.L. [7]), and the corresponding Higgs-mediated contribution to the previous LFV processes 
can be relevant, even for large soft SUSY masses. 

In the previous related literature there are, to our knowledge, just a few theoretical 
computations of some of these LFV semileptonic r decays induced by SUSY loops. In 
particular, r — > /irj was first computed in [31] within the context of the unconstrained MSSM 
and in the approximation of large tan/5. A more refined analysis of this channel, r — > /irj', 
r — > /!7r, and t — > jip was done in [32] for the unconstrained MSSM scenario and large 
tan/3 approximation as well, but they used an effective lagrangian framework for the LFV 
operators. An estimate of r — > /irj with the use of the mass insertion (MI) approximation 
for the relevant lepton flavour mixing parameter between the r and /x sectors, 632, has been 
performed in [33]. The decay mode r — > fiK + K~ has been estimated in [34] within the mass 
insertion and leading logarithmic (LLog) approximations for £32, and taking into account 
only the Higgs-mediated contribution in the large tan/5 limit. In all these previous works no 
connection with the neutrino sector was considered and the hadronisation of quark bilinears 
in the final state is simply parameterised in terms of the proper meson decay constants and 
meson masses. Other estimates of some of these LFV semileptonic r decays in different 
contexts, like >SO(10)-SUSY-GUT model with universal soft breaking terms [35] and Littlest 
Higgs model [36], have also been performed in the literature. 

Our analysis presented here is more complete in several aspects. First, we include both 
Z-boson and A°-boson mediated contributions to r — > jiP (P = i],i]',7r ), and both 7 and 
H°, /i°-bosons mediated contributions to r — > fiK + K~ . The other channels, r — > {iK°K° and 
t — > fi^n have not been estimated previously. We include 7 and H°, /i°-bosons mediated 
contributions in r — > fiK°K°. The case r — > fi7T 7r° can only be mediated by H°, h°- 
bosons. Second, we do not use either the mass insertion nor the LLog approximation and 
our analytical computation is valid for all tan/5 values. Third, we make a connection with 
neutrino physics by requiring compatibility through all this work with the neutrino data for 
masses and mixing angles. Fourth, we perform the hadronisation of quark bilinears with close 
attention to the chiral constraints, guided by the resonance chiral theory [37] that has proven 
to be a robust framework for the analyses of hadrodynamics when resonances are involved. 
The 7 amplitude, due to its pole at q 2 = 0, is most sensitive to the hadronisation procedure. 
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Hence the hadronisation of the electromagnetic current, that drives the 7 contributions, has 
been carried out by a careful construction of the vector form factor that matches both the 
chiral low-energy limit and the asymptotic smoothing at high q 2 [38]. Those final states 
driven by heavy intermediate bosons like the Z° or Higgses, on the other side, do not require 
such an involved scheme. In these cases we have used the leading chiral approximation of 
Chiral Perturbation Theory that we know, for sure, it has to be fulfilled by the hadronisation. 
The advantage of our approach is that it provides the most successful description up to date 
of the hadronic tau decays and it can be systematically improved by further developments 
of the appropriate form factors, whether axial-vector, scalar or pseudoscalar cases. 

The rest of this paper is organised as follows. The theoretical framework for the compu- 
tation of LFV semileptonic r decays is described in Section 2. This includes a short review 
of the SUSY-seesaw scenarios that we work within, CMSSM and NUHM, and a description 
of our procedure for hadronisation of quark bilinears within the context of Chiral Pertur- 
bation Theory (%PT) and Resonance Chiral Theory (R%T). In Section 3, the analytical 
results of the full one- loop branching ratios BR(r — > fiPP), BR(r — > /iP), BR(r — > fip) and 
BR(r — > are presented. Section 4 is devoted to the numerical results and discussion. It 
includes, in addition a comparison between the full one-loop and approximate results. A set 
of useful approximate formulas for the semileptonic tau decay rates that are valid at large 
tan/5 are derived. A critical comparison with previous predictions in the literature is also 
included in this Section 4. Finally, Section 5 summarises the conclusions. 



LFV semilep. r decays 



BABAR Belle 
1.5 x 1(T 7 [27] 6.5 x 1(T 8 
1.4 x 1(T 7 [27] 1.3 x 1(T 7 [26] 
1.1 x 10~ 7 [27] 1.2 x 1(T 7 [26] 
2.0 x 10~ 7 [25] 
1.3 x 10~ 7 [39] 
4.8 x 10~ 7 [25] 



BABAR & Belle 



5.1 x lO" 8 [15] 
5.3 x 10~ 8 [15] 
5.8 x 10~ 8 [15] 



BR(r -> fir]) 
BR(r -> fir]') 
BR(r -> fm) 
BR(r -> fip) 
BR(r -> 11(f)) 
BR(r -> fin+n-) 
BR(r -> ^ttV) 
BR(r -> pK + K-) 
BR(r -> jj,K°K°) 



8.0 x 10 



-7 



Tabic 1: Present upper bounds for LFV semileptonic r decays. 
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2 Framework for LFV semileptonic r decays 



In this section we describe the theoretical framework for the computation of the LFV semilep- 
tonic r decay rates. First we present the scenario for the generation of LFV in the r-/i sector, 
then we summarise the main ingredients to perform the hadronisation of quark bilinears 
within the context of yPT and RyT. 

2.1 LFV in the SUSY-seesaw scenario 

The SUSY-seesaw scenario that we work within contains the full MSSM spectra and, in 
addition, the three right handed neutrinos and their SUSY partners. It is defined in terms 
of both the SUSY and neutrino sector parameters which are summarised in the following. 

Regarding the SUSY sector we choose to work in two different constrained MSSM scenar- 
ios: The usual Constrained MSSM (CMSSM) with similar input parameters as in mSUGRA 
models, and the so-called Non-Universal Higgs Mass (NUHM) scenarios with two additional 
parameters defining the non-universal soft Higgs masses. The corresponding sets of input 
parameters in these two scenarios are: 



where M , M 1 / 2 and A are the universal soft SUSY breaking scalar masses, gaugino masses 
and trilinear couplings at the gauge coupling unification scale, M x = 2 x 10 16 GeV. Notice 
that M and A define also the soft parameters in the sneutrino sector. The other CMSSM 
parameters are, as usual, the ratio of the two Higgs vacuum expectation values, tan (5 = i^/fi, 
and the sign of the ji parameter, sign(/i). The departure from universality in the soft Higgs 
masses of the NUHM is parameterised in terms of two parameters 5\ and <5 2 by, 



Notice that with the choice <5i i2 = one recovers the universal case defined by the CMSSM 
scenario. For simplicity, and to further reduce the number of input parameters, in all the 
numerical estimates of this work we will take M = My 2 = Mjusy, A) — and sign(/x) = +1. 

To evaluate these two sets of parameters at low energies (taken here as the Z gauge boson 
mass m z ) we solve the full one-loop Renormalisation Group Equations (RGEs) including the 
extended neutrino and sneutrino sectors. For this and the computation of the full spectra 
at the low energy we use here the public FORTRAN code SPheno [40] . 



CMSSM : M , Mi/a , A) , tan (5 , sign(/x) , 
NUHM : M , M 1/2 , A , tan (3 , sign(/x) , M Hl , M H2 , 



(1) 



M 2 Hi = M^(l + 5 1 ), M| 2 = M 2 (l + 5 2 ). 



(2) 
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Regarding the neutrino sector, we use the seesaw mechanism for neutrino mass generation 
which is implemented here to the case of three right handed neutrinos. The usual input 
parameters in this case are the three right handed Majorana masses, Mr 123 and the neutrino 
Yukawa coupling 3x3 matrix, Y v . The Dirac neutrino mass matrix is then related to the 
Yukawa couplings by m D = Y v v 2 , where v 2 — v sin/3 and v = 174 GeV. In this seesaw 
scenario, the physical Majorana neutrinos consist of three light ones, ^1,2,3, with predicted 
masses being typically m Vl23 ~ 0(m 2 D /Mji), and three heavy ones, ^2,3, with masses 
mjvi 2 ,3 ~ Mr 12 3 . However, instead of this we will use another set of input parameters which 
are more convenient to accommodate the experimental data on light neutrino masses and 
generational mixing angles. Within this parameterisation, the Dirac mass matrix and, hence, 
the Yukawa coupling matrix, are derived in terms of the physical neutrino masses, neutrino 
mixings and a generic complex orthogal 3x3 matrix, R, as follows [41], 

m D =Y u v 2 = iyj m% ag R \j m^ iag U MNS , (3) 

where, 

m^ ag = diag (m Nl , m N2 , m N3 ) , (4) 

m^ iag = diag {m Ul ,m U2 ,m V3 ) , (5) 

and we use the standard parameterisation for the unitary matrix Umns [42, 43] containing 
the three generational mixing angles 9u, #13 and #23 and the three CP violating phases, 5, 
0i 5 2- In turn, the R matrix is parameterised in terms of three complex angles, 9i (i — 1,2, 3) 
as [41] 

< C 2 C 3 -Ci S3 - Si S 2 C 3 Si S3 - Ci s 2 c 3 \ 

R = c 2 s 3 ci c 3 - si s 2 s 3 -si c 3 - ci s 2 s 3 , (6) 

\ S 2 SiC 2 CiC 2 J 

with q = cos#j and Sj = sin^j. One interesting aspect of this matrix is that it encodes 
the possible extra neutrino mixings (associated with the right-handed sector) in addition to 
the ones in £/mns- Notice also that the previous Eq. (3) is established at the right-handed 
neutrino mass scale Mr, so that the quantities appearing in it are indeed the renormalised 
ones, namely, m^ iag (Mr) and ?7mns (Mr). These latter are obtained here by means of the 
RGEs and by starting the running from their corresponding renormalised values at mz, 
m dia g ( mz ) anc i C/mns ( m z) which are identified respectively with the physical m^ iag and 
Umns from neutrino data. 

Concerning our choice for the size of the physical neutrino parameters, we shall focus in 
this work on scenarios where both light and heavy neutrinos are hierarchical, 

m Ul < m V2 < m V3 , 
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m Nl < m N , 2 < 1TIN3 , 



and set the numerical values for the light neutrino parameters to the following ones which 
are compatible with present data [7, 44] 

m\ ~ , m 2 2 = Am s 2 ol = 8 x 1(T 5 eV 2 , m\ = Am 2 tm = 2.5 x 1(T 3 eV 2 , 

012 = 30° , 23 = 45° , 013 ^ 0, 5 = 0! = 2 = . (7) 

Notice that, for simplicity, the three CP violating phases, have been set to zero. We have 
also set to zero the 0i3 mixing angle and the lightest neutrino mass in order to minimise as 
much as possible the LFV in the \i — e sector. In fact, we have checked that for the explored 
parameters region in this work, this \i — e LFV is below the sensitivity of the present data 
from ji — > cy, /i — » 3e and fj, — e conversion in nuclei. 

In summary, the input parameters of the neutrino sector for the present work are: 

Seesaw : rriN 12 _ 3 , 0i,2,3 • (8) 

Regarding the generation of LFV in these constrained MSSM scenarios, we remind that 
all lepton flavour mixing originates solely from the neutrino Yukawa couplings. These Y u 
first induce flavour violation in the slepton sector by the RGE running of the soft SUSY 
breaking parameters from Mx down to the electroweak scale mz- It is manifested in the 
non- vanishing values of the off-diagonal elements of the slepton squared mass matrix at mz- 
We perform this running by solving the full set of one-loop RGEs including the neutrino 
and sneutrino sectors. The resulting slepton mass matrices at mz are then diagonalised and 
the previous flavour mixing is then transmitted to the mass eigen-values and eigen-states. 
Therefore, in this work where we deal with physical states, all flavour mixing is implicit in the 
resulting physical charged slepton masses, m~ , ■■,m'~ , sneutrino masses, m% l , fnl 2 , m~ 3 and 
the corresponding matrices that rotate from the electroweak to the slepton and sneutrino 
mass bases, respectively, it!' and R v . The LFV in the physical processes, like lj — > 1^, 
lj — > 3/j, n — e conversion in nuclei, and the semileptonic r decays studied here, are then 
generated by the SUSY one-loop contributing diagrams that contain these slepton physical 
masses in the internal propagators, and also the previous rotation matrices in the interaction 
vertices, which connect between different lepton generations. A complete set of Feynman 
rules for the relevant LFV vertices can be found in [21,24]. 

Finally, in order to illustrate more quantitatively how important can be the size of the 
flavour mixing between the stau and smuon sectors, in the CMSSM-seesaw scenario, we 
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include next the predictions of the mixing parameter £32 that is defined in the LLog approx- 
imation as, 

1 (3M 2 + A 2 ) 



'32 



07T ^ W SUSY 



W^) 32 , (9) 



where L is a 3 x 3 diagonal matrix whose elements are, La = Xog^Mx/mN^ and Msusy is an 
average SUSY mass. This phenomenological parameter <5 32 measures the amount of flavour 
mixing between the second and third slepton generations in the left-handed sector (LL), 
which is by far the dominant one. The corresponding mixing in the right-handed slepton 
sector is extremely suppressed by the smallness of the lepton masses which appear as global 
factors in the definitions of those (RR and RL) mixings (see, for instance, [21]). 

One can estimate £32 from the previous parameterisation of the seesaw model in Eq. (3) 
by simply plugging in Eq. (9) the value of (Yj L Y v ) 32 from the following expression, 

V 2 ( Y u LY ^) 32 = ^33^3 [(v / ^ C l C 2 C 13C23 - v 7 ^ 5 * ^^^) 

( v /m^"cic 2 s 2 3 + ^frn^s ^20^023)] 
+ L 22 m N2 [{^/m^{-s\cl - cls* 2 s* 3 )c 2 3 + ^/m^ 2 (s* v s* 2 sl - CiCg)^^) 

( v /m^"(-SiC 3 - C1S 2 S 3 )S23 + y/m^(ciC 3 ~ S^^)^)] 

+ Lu m Nl [(y^m^islsl - c\s* 2 cl)c 2 3 + ^/m^{.sls* 2 cl + clsl)c 12 s 23 ) 

(v / ^( S l S 3 - C1S 2 C 3 )S23 - V / ^"( S l S 2C 3 + CiS 3 )Ci 2 C2 3 )] , (10) 

where, Sij = sin^- and c^- = cos%, and we have already set m Vl = and #13 = 0. 

The numerical predictions for |Yj, 32 | and |<5 32 | as a function of 9 2 are shown in Fig. 1. Here 
we have set M SUS y = M = M 1/2 = 250 GeV, A = 0, (m Nl ,m Na ,m N3 ) = (10 10 , 10 11 , 10 14 ) 
GeV, and the light neutrino parameters are those in Eq. (7). We see clearly that follows 
the same pattern as \Y^ 2 \ (and (Yj 53 !) and can reach large values in the range 0.1-1 for several 
choices of \9 2 \ and arg(6> 2 ). Notice also that the predictions for |5 32 | corresponding to Yukawa 
couplings larger than about 4 are not shown, because through all this work perturbativity 
in all the gauge and Yukawa couplings are imposed. This is set numerically in the Spheno 
program by the requirement |Y iy | 2 /(47r) < 1.5 and corresponds to a maximal predicted value 
of about \5 32 \ < 0.4. The corresponding predictions with respect to 9\ are very similar to 
those of 92 and are not shown for brevity. The value of {632 | is practically independent on 
9 3 . For the rest of this work we will set #i i3 = and use just 9 2 as input parameter. 

The numerical predictions for | ^32 1 as a function of the heaviest neutrino mass, m^ 3 
are shown in Fig.2. |5 32 | values within the range 0.1-1 are obtained for large values, 
say within the interval 10 13 — 10 15 GeV. Notice that the predictions enter into the above 



8 



*********** 
± "+^+++++++^+++++- 



f++++ 



+++++ +H : 



M = M 1/2 = 250 GeV, A = 
tan p = 50, 8, = e 3 = 
5, = 8 2 = 



0.5 



1.5 



|6 2 | 



arg e 2 = • 

arg 8 2 = ji/8 + 

arg 8 2 = jt/4 * 

arge 2 = 3n/8 a 

arg B 2 = it/2 o 



2.5 



10 u 



10 1 



10"' 



10" 1 



10' 



10' 



*»»«»88Sllll********* 

++ +++++. 



***"* 
++++++++++ ++ " 



+ ++++ ++ 



+++ 



+++■ 



m N = (10 10 , 10 11 , 10 14 )GeV • 
M = M 1/2 = 250 GeV, A = 
tan p = 50, 8, = e 3 = 
8, = 8 2 = 



0.5 



1.5 



|e 2 | 



.++■ 



arg 6 2 = ■ 

arg 6 2 = it/8 + 

arg 6 2 = jc/4 * 

arg 6 2 = 3it/8 a 

arg 6 2 = nl2 o 



2.5 



Figure 1: \Y^ 2 \ and [632], in the CMSSM-seesaw scenario, as a function of \02\, for arg#2 = 
{0, 7r/8, tt/4, 37r/8, 7r/2} (dots, crosses, asterisks, triangles and circles, respectively). Both \02\ and argf?2 
are given in radians. The predictions for |Y„ 33 | are practically indistinguishable from those for |Y^ 2 |. 
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Figure 2: |5 32 |, in the CMSSM-seesaw scenario, as a function of mjv 3 



commented non-perturbative region for values larger than m^ 3 = 10 GeV, and for the 
particular choice of O2 = 2.9 e^. Concretely, the value | ^32 1 = 1, which is interesting for 
later discussion and comparison with other works, corresponds to m^ 3 = 3 x 10 14 GeV and 
lies clearly in the non-perturbative region. Finally, just to mention that is not much 
dependent on tan/? nor onmjv 12 . The values of these two heavy neutrino masses will be set 
in the following to the reference values 2 = 10 10 , 10 11 GeV. 
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2.2 Hadronisation of quark bilinear currents 

Semileptonic decays of the tau lepton are a relatively clean scenario from the strong interac- 
tion point of view. Hadrons in the final state stem from the hadronisation of quark bilinears, 
namely ip Tip, where ip is a vector in the SU (3) p flavour space and Y is, in general, a matrix 
both in the spinor and the flavour space. 

An appropriate framework to handle the procedure of hadronisation is provided by the 
large- N c expansion of SU(N C ) QCD [45], being N c the number of colours. In short it 
stays that in the Nc — > oo limit any Green function is given by meromorphic expressions 
provided by the tree level diagrams of a Lagrangian theory with an infinite spectrum of 
zero- width states. Though we do not know how to implement fully this limit, a fruitful [46] 
if debatable [47] approach lies in cutting the spectrum, keeping only the lightest multiplets 
of resonances. We will attach to this tenet as a guiding principle. 

A suitable tool to realise the 1/N C expansion is provided by chiral Lagrangians. In 
those processes where hadron resonances do not play a dynamical role, %PT [48,49] is the 
appropriate scheme to describe the strong interaction of Goldstone bosons (tt, K and rj). 
This is the case, for instance, of r — > \iP (being P short for a pseudoscalar meson). When 
resonances participate in the dynamics of the process, as in r — > fiPP, it is necessary to 
include them as active degrees of freedom into the Lagrangian as it is properly done in the 
R/\T frame [37]. Hence we will make use of R/\T, that naturally includes x?T, to hadronise 
the relevant currents that appear in the processes under study here. 

We consider bilinear light quark operators coupled to external sources and added to the 
massles QCD Lagrangian : 

£ Q cd = £q C d + q [7m K + 75 a M ) - ( s - ipis)] q , (11) 

where vector (v^ = vf\ l /2), axial- vector (a M = af A l /2), scalar (s = SjA*) and pseudoscalar 
(p = Pi\ l ) fields are matrices in the flavour space, and £q CD is the massless QCD La- 
grangian 1 . This Lagrangian density gives the QCD generating functional ^qcd [v,a, s,p\ 
as 

e iZ QC D[v,a,s,p] _ J [JJG^Jjq^Jjqj^ J d4xC QCo[<l,q,G,v,a,s,p] ^ ^2) 

In order to construct the corresponding Lagrangian theory in terms of the lightest hadron 
1 The Gell-Mann matrices A* are normalised as (A,Aj) = 25y and the gluons are denoted here by G M . 
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modes we need to specify them. The lightest U (3) nonet of pseudoscalar mesons 

4>{x) 



a=0 
/ 1 



(13) 



V2 



7T° + 



V 



1 

K~ 



1 



71 



1 



"7T + 7^ 8 + 7?° 

K° 



K+ 



\ 



is realised nonlinearly into the unitary matrix in the flavour space : 

[ $ 

u((fi) = exp i -j^ 

Hence the leading 0(p 2 ) %PT SU(3) L ® SU(3) R chiral Lagrangian is 



(14) 



(15) 



where 



X+ = u^xu^ + ux^u , x = 2B o(s + ip) , 



(16) 



and (. . .) is short for a trace in the flavour space. Interactions with electroweak bosons can 
be accommodated through the vector = (r^ + i^)j2 and axial- vector a M = (r^ — £p)/2 
external fields. The scalar field s incorporates explicit chiral symmetry breaking through the 
quark masses s — M. + ... and, finally F ~ F n ~ 92.4 MeV is the pion decay constant and 
BqF 2 = — (0\ipifj\0) o in the chiral limit. The chiral tensor x provides masses to the Goldstone 
bosons through the external scalar field, as can be seen in Eq. (16). Indeed in the isospin 
limit we have : 



X = 2B M + ... = 
Hence we identify : 

B m u 



mi 



2m\ 



m 



\ 



2 / 

IT / 



+ 



(17) 



B m d 



1 2 



2 Notice that though we include a U (3) nonet we are not relying on the U (3)l®U(3)r chiral Lagrangian [50] 
on grounds of predictability, as the latter introduces unknown functions. 
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B m s = m 2 K -^m 2 n , (18) 

that will be useful when considering the Higgs contributions. The mass eigenstates rj and rj' 
are defined from the octet r)$ and singlet rj states through the rotation : 



cos — sin 9 
rj' I \ sin 9 cos 9 




(19) 



and we input 3 a value of 9 ~ — 18°. 

The hadronisation of a final state of two pseudoscalars is driven by vector and scalar 
resonances though the latter, because their higher masses, play a lesser role and we will not 
include them in the following. We will introduce the vector resonances in the antisymmetric 
formalism; hence the nonet of resonance fields [37] is defined by analogy with Eq. (13) 
with the same flavour structure. By demanding the chiral symmetry invariance the resonance 
Lagrangian reads : 

r _ rV i t 

'(2) , 

where 



C v = £L + CY 2) , (20) 



1 M 2 

4) = ^^D + ^^O . (2i) 

and in the latter the subscript (2) indicates the chiral order of the tensor accompanying V^ v . 
In Eq. (21) we have used the definitions : 

V,X = d,X+[T^,X} , (22) 
r M = -[u^dp-irju + ufip-ie^v!] , 
= u F£'u* + u*Fg'u, 

being F^ U R the field strength tensors associated with the external right and left fields. The 
couplings F v and Gy are real. 

Accordingly our R%T framework is provided by : 

£ RxT = £f + C v , (23) 



3 The values of 9 in the literature range between 6 <~ —12° up to 9 ~ —20° [51]. 
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and the contribution of the low modes to the QCD functional is formally given by : 



J Z QC v[v,a,s,p] _ J [£) u ] [DV] e* / d4:r£ RxT[«.V,v,a,s,p] ^4) 

low modes 

With this identification we can already carry out the hadronisation of the bilinear quark 
currents included in Eq. (11) by taking the appropriate partial derivatives, with respect to 
the external auxiliary fields, of the functional action, 



S* = -qXq 



dsi 



3=0 



3=0 



Ai - ^ 9C RxT 



3=0 



P i = q ^75 A* q 



dpi 



(25) 



3=0 



where j = indicates that all external currents are set to zero. This gives : 



V 1 = 
V 

S l = 

p i = 



4 

F 2 



( X (u u' — v) uj ) - 
( X 1 (uu^v) + v) u^u)) , 



2^/2 



( X d u {u^V^u + uV v ^)), 



-B F 2 (X i (uV + uu)) , 
- B F 2 ( X (uV — wiij ) . 



(26) 



With these expressions we are able to hadronise the final states in r — > fiPP and r — > fiP 
processes as we explain now : 

7 contribution 

The photon contribution to the decay into two pseudoscalar mesons is driven by the elec- 
tromagnetic current : 

u,d,s 



(27) 



where Q q is the electric charge of the q quark in units of the positron charge e. The electro- 
magnetic form factor is then defined as : 

(Pi(pi)P 2 (p 2 )|V; em |0) = (pi- P 2)»F? P2 (s), (28) 

where Fy lP2 (s) is steered by both 1 = 1 and 1 = vector resonances, in particular the 
p(770) that is the lightest of resonances. Due to the q 2 = pole of the photon propagator 
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this is, by far, the dominant contribution to this hadronic final state. Hence the result is 
more sensitive to the construction of this form factor. Accordingly we will elaborate a more 
complete expression than the one provided by the vector current in Eq. (26) though it will 
reduce to this one in the Nc — > oo limit, including only one multiplet of resonances and at 
q 2 <C M 2 A proper construction of Fy lP2 (s) is given in Appendix B. 

Z° contribution 



Here both vector and axial-vector currents do contribute. In terms of the quark fields these 
are : 



jZ yZ , \Z 

,rZ 9 - 

^ 2 cos W 



2 sin 2 9 W Q - T 3 (9) 



q 



with Q = diag(2, — 1, — 1)/3 and = diag(l, — 1, — 1)/2 the electric charge and weak 
hypercharges, respectively, g is the SU (2) gauge coupling and 9 W is the weak angle. 

In order to proceed to the hadronisation of these currents one has to write the currents 
in Eq. (29) in terms of V* and A 1 ^ defined in Eq. (26). This gives 

2 sin 2 ^ (Q {uu^v) — u^u^u)) — (T^ {uu^v) — u^u^uj) 



9 F 2 



2 cos 9w 2 

?2 



K = TT-Z-n- 4- ( T 3 (9) (tttt„«t + u ^ u ) ) . (30) 
L COS v\\r Z 

Notice that the vector current contributes to an even number of pseudoscalar mesons while 
the axial-vector current provides 1,3,... mesons. 

Higgs bosons contribution 

Hadronisation of scalar Higgs bosons like h° and H° into two pseudoscalar mesons proceeds 
through the scalar current while the pseudoscalar A Higgs boson hadronises through the 
pseudoscalar current into one pseudoscalar meson. As Higgses are rather massive the hadro- 
nisation is not so sensitive to resonances as in the case of the photon contribution. Hence 
we will not elaborate on scalar of pseudoscalar form factors (analogous to the vector case 
defined by Eq. (28)) that, moreover, are not so well known. We will rely in the following 
scalar and pseudoscalar currents, 

uTu = ±J 3 + -Lj 8 + -Lj°, 
2 2V3 V6 
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Figure 3: Contributions to LFV semilcptonic r decays into one and two pseudoscalar mesons. 



dT d 



sY s 




(31) 



where T = — 1 for J 1 = S 1 *, V = 275 for J % = P % and the S l and P % currents are given in 



3 Analytical results of the LFV semileptonic r decays 

In this section we present the analytical results of the branching ratios for the LFV semilep- 
tonic r decays: r -> fiPP, with PP = tt+tt" ,7r°7r°, K + K~, K°K° and r -> \iP y with P = 
it, i] and v( . The predictions for the r — > /xp and r — * fup channels, which are related to 
r — > fi7T + 7r~ and r — > fiK + K~ , fiK°K° respectively, will also be included. 

3.1 Predictions for r — ► //PP 

The semileptonic r — > pPP channels can be mediated by a photon, a Z gauge boson and a 
CP even Higgs boson, h° and if . The various contributing diagrams are depicted in Fig. 3. 



Eq. (26). 
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In these diagrams, the LFV vertex is represented by a black circle and the hadronic vertex 
by a white box. The Z-mediated contribution is expected to be much smaller than the 7- 
mediated contribution due to the 0{l/m 2 z ) suppression factor in the amplitude from the Z 
propagator. This has been shown to happen in the leptonic channels like r — > 3/i, where the 
Z-mediated contribution to its branching ratio has been estimated to be a factor 1CT 3 — 1CT 5 
smaller than the 7-mediated contribution, for tan/3 = 5 — 50 [21]. Consequently, we have 
neglected here the Z contribution to the r — > \xPP decays. By using again this comparison 
with r — > 3/i, the 7 contribution to r — > jiPP is expected to be the dominant one, and the h° 
and if °-mediated contributions are expected to be relevant only at large tan (3. Therefore, 
we have included these three 7, h° and H° contributions in the computation. 

The total amplitude for the r — > u,PP process can be written as, 



T = 



H 



(32) 



where T 7 and Th = T h o + T H o are the amplitudes of the 7-mediated and if -mediated con- 
tributions respectively. First we present the result of T 7 and Th in terms of the final state 
quarks, that is for r — > fjiqq, and in terms of the corresponding r-/i LFV form factors: 

Jl [k 2 lfl [A\P L + AfP R ) + im T a^k u {A L 2 P L + A*P R )} r 
e 2 Q q _ 



x 



k 2 



q QYQ : 



(33) 



T H = 



h°,H° h p 

+ H^S^ q [jlP L r] [qP R q] + H%>S® [JIP R t] [qP L q]} , 



(34) 



where, k is the photon momentum, Q q the electric charge of the quark q in units of the 
positron charge e, Pl,r = (1 -F 7s)/2, ^ t is the r lepton mass, and m^o, m#o are the Higgs 
boson masses. Notice that the momentum of the Higgs propagators has been neglected 
against the Higgs boson mass. The Higgs boson couplings to quarks are correspondingly 
given by, 



S 



(p) 

Lm 



9 



-a 



(p)* ' 



2m 



w 



sin (3 



m„ 



c(p) 
°L,(d,s) 



9 



(7i 



2m 



w 



cos (3 



m d . 



c(p) 



q(p)* 



(35) 



where m g is the q quark mass, the W gauge boson mass, g the SU (2) gauge coupling, 
and 



» 




a ip) 



f 
\ 



COS (X 

sin a 

-icos/9 J 



(36) 
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The three entries for the index (p) in the previous expressions and in the following ones cor- 
respond to Hp = h°, H°, A , respectively. The angle a rotates, as usual, from the electroweak 
neutral Higgs basis to the mass eigenstate basis. 

The LFV form factors in Eq. (33) describe the effective 777/ vertex and get contri- 
butions from the SUSY one- loop diagrams depicted in Fig. A.l. The full results for these 
form factors can be found in the literature [2,21] and are collected in Appendix A.l for 
completeness. Notice that we are presenting all the results in the physical mass eigenstate 
basis for all the particles involved. Therefore the LFV is encoded in the physical slepton 
and sneutrino masses and in the corresponding slepton and sneutrino rotation matrices. The 
later appear in the chargino-sneutrino-lepton and neutralino-slepton-lepton couplings. Sim- 
ilarly, the LFV form factors H^ R in Eq. (34) describe the effective H p rfi vertex and get 
contributions from the SUSY one-loop diagrams shown in Fig. A. 3. These set of diagrams 
where computed in [52] and the results are collected in Appendix A. 3. Again the LFV is 
encoded in the slepton and sneutrino masses and in the rotation matrices. 

The next step is to hadronise the quark bilinears appearing in Eqs. (33) and (34). For this, 
we proceed as explained in Sec. 2.2. The quark bilinears in T 7 , [57^?], are hadronised via the 
electromagnetic current in Eq. (27) which, for the final state with two mesons Pi(pi)P2(p2), 
is then written in terms of the corresponding electromagnetic form factor, Fy lP2 , by means 
of Eq. (28). Thus, one gets the photon amplitude in terms of the final state hadrons: 

T ~< = TS F v P '( e )~P [* 2 *i-« K p t + A ? p R) + 2rm T K<VP2 (^Pl + A*P R )] t. 

(37) 

The expressions of the F PlP2 form factors for each of the final states, P1P2 =7r + 7r~, K + K~ 
and K°K° are collected in Appendix B. Obviously, the 7r°7r° final state does not get photon- 
mediated contributions since 7 does not couple to 7r°7r°. Hence we set Fy n ° = 0. 

The quark bilinears in Th, [qPl,rq\, when hadronised in a final state of two mesons, get 
contributions just from scalar currents, S\ but not from pseudoscalar currents, P % . Then, 
one substitutes \<}Pl,rq\ by [(— 1/2)(— qq)], where (—qq) is given in Eq. (31), and the relevant 
scalar currents, S°, S 3 and S* 8 , are written in terms of two mesons by using Eq. (26). This 
gives: 



S 3 = -£ ( 

cs _ B o 
s ~ 7=3 



— (cos 9- V2sin0) A + -y= (v^costf + sinfl) n rj' + K + K~ - K°K° 
K + K- + K°K° - 27T+7T- - ttV + (cos 2 6 + 2^2 sin d cos rjr] 
+ 2 (v^sin^ + sinflcosfl- v^cos 2 ^ 7777'] , 
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S° = -Bo^ [2n + 7i- + 2K+K- + 2K°K° + ttV + rjrj] . 

Thus, one gets the Higgs boson amplitude in terms of the final state hadrons: 



(38) 



Th = Yl V 

p=h°,H° 



(p) 



T, 



(39) 



where 



^pp 



9 



2mw 2mj I 



9 



d {p) = 
pp 2m w 2m\ 



l — (j£\pp) + j$\ppj) (hV + hP) , 

Hp 

L- (j«(PP) + jW(PP)) (lf« - 1#>) , 



(40) 



and 



jf(^K-) 

4 p) (^°i?°) 



._(p)* ' 

°2 

sin /3 



m 2 + 



(p)* 



cos /3 



1 

4 

1 / a 



—a 



(p)* ' 



sin /3 



(p)* 



cos /5 



(41) 



2 \ cos /3 

4 p) (pp) = 4 p) *(pp). 

Notice that in Eq. (41) we have already used the relations between the quark and the meson 
masses of %PT given in Eq. (18). 

Finally, we get the following result for the branching ratio: 

!~ f Smax /'Wax 1 

BR ^" PP > = 64^/ " S I dt 2^ T ^ < 42 » 

i ' ' J d min ^min A ¥ 

where T T is the total r decay width, and the coefficient Kpp is 1 for PP = 7r + 7r~, K + K~ , K°K° 
and 1/2 for PP = 7r°7r°. In addition 

1 r, 



4s 



(ml -ml) 2 - (A 1/2 (s, mp, m 2 p) =F A 1 / 2 (m 2 ,s,m 2 ))' 



s min = Amp , s max = (m T - m^f , A(x, y,z) = (x + y- zf - 4xy . 
The averaged squared amplitude is, 



1 

8 m, 



^i(s) + g 2 (s)t + 9 3 (s)t 2 ] 



(43) 



(44) 
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where 

9i{s) = h + h x s + h 2 s 2 + h 3 s 3 , 

5-2(5) = jis + j 2 s 2 + j 3 s 3 , 

g 3 (s) = hs + k 2 s 2 , (45) 

with 

ho = -8M 2 P m 2 T (ml-m 2 T ) 2 (A 2 A 2 * + A+Ap) 

+ 2 (m M + m T f c H c* H + 2 (m M - m T ) 2 , 
h x = -8m 2 . (m r m M + M 2 pf A 2 A 2 * - 8m 2 T (m r m M - Mp) 2 A+A+* 

+8 (m M - m T ) m T (m M + m T ) 2 Ml (A±*A 2 + AfA^*) 

+8 (m M - m T ) 2 m r (m M + m T ) Mj, (A|*A^ + A^Ap) 

-2 (m M + m T ) (m 2 + m 2 T + 2M 2 P ) (c H A+* + c* H A+) 

+2m T (ml + m 2 T + 2M 2 P ) (c H Ap + c* H A+) 

-2 (m M - m T ) (mj + m 2 + 2M 2 P ) (d H A^* + rf^Af) 

+2m r (ml + m 2 T + 2M 2 P ) (d H A 2 * + d* H A~ ) 

-2c H c* H - 2d H d* H , 
h 2 = 2 \(ml + m 2 T ) 2 + 4M P + 8m fl m T M P ] A+A+* 

+2 Uml + m 2 ) 2 + 4MJ, - 8m M m r M|] A[A[* 

+ 2m 2 [(m M - m T ) 2 + 4M|] A\A\* + 2m 2 [(m M + m r f + 4M|] A 2 A 2 * 

- 2m T (m M - m r ) [(m„ + m T f + 4M|] + A~A~ *) 

- 2m r (m„ + m T ) [(m M - m T f + 4M|] + A\Ap) 

+ 2 (m M + m T ) (c H A+* + c* H A+) + 2 (m M - m r ) (rf^Af* + d* H A{) 

- 2m T [c H Aj* + c* H A+ + d^A-* + cf^j] , 

h 3 = -2 (m M - m T ) 2 A^Ai* - 2 (m M + m r ) 2 A+A+* - 2m 2 [AT A 2 * + A+A+*] 

+ 2m r (m M - m T ) [A^A 2 * + *A 2 ] + 2m T (m M + m T ) [A+A+* + A+*A+] , 
3l = 8m 2 (ml + m 2 +2M 2 P )(A 2 A 2 * + A+A+*) 

- Am T [c H A$* + c* H A+ + 44* + d* H A 2 ] 

+ 4 (m M + m T ) (c#A+* + c* H At) + 4 (m M - m T ) (rf^Af * + d* H A^) , 
h = -8 (ml + m 2 T + 2M 2 P ) [A+A+* + A^A^] - 8m 2 T [A+Ap + A 2 A 2 *} , 
j 3 = 8 (A-.A-p + , 
fci = -8m 2 (A 2 A 2 * + A+ A+*) , 
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T A o = 



h = 8 (A^r* + A+At*) , (46) 

and 

Af = f s F' p (s) {A? ± A?) , c H = cfj + cfP , d H = dg? + cff . (47) 
3.2 Predictions for r — > //P 

The semileptonic r — > /iP channel can be mediated by a Z gauge boson and a CP odd 
A Higgs boson, as represented in Fig. 3. Both contributions are included here. The total 
amplitude for this r — > \xP decay can then be written as, 

T = T z + T A o , (48) 

where Tz and T40 are the Z and A mediated amplitudes respectively. As in the previous 
case, these are first evaluated in terms of the final state quarks, that is for r — > fxqq, and in 
terms of the corresponding t — fi LFV form factors: 

Tz = \jlh,{F L P L + F R P R )]r .q\r(z^P L + Z^P R )]q, (49) 

J- K 0) < } ^r] [?P L ?] + HfW£ [HP R r] [qPm] 

+ HfW£> [PP L t] [qP R q] + H^S^ [pP R r] [qP L q]} , (50) 

where Z^ = (—g/ cos6\y)(T^ — Q^sin 2 ^) and Z^ — (g/ cos9w)Q q sin 2 6\y are the Z 
couplings to quarks, and ^ and S R A "* are the A couplings to quarks, which are given by 
the third entry in Eqs. (35) and (36). Notice that, as in the previous cases of h° and H°, we 
have neglected the k 2 in the Z and A propagators. 

The LFV form factors Fl )R in Eq. (49) describe the effective Zt/j, vertex and receive 
contributions from the SUSY one-loop diagrams depicted in Fig. A. 2. The results for these 
form factors where found in [2] and corrected in [21]. We collect them in Appendix A. 2, for 
completeness. The LFV form factors in Eq. (50) describe the effective A°Tfi vertex and, 
as in the previous Hrfi vertices with H = h°, H°, receive contributions from the one-loop 
diagrams in Fig. A. 3. The corresponding results are collected in Appendix A. 3. 

The hadronisation of the quark bilinears in Tz proceeds by means of the vector and axial- 
vector currents in Eq. (29), which in turn are written in terms of one P meson by means of 
Eq. (30). This leads to: 

V* = 0, (51) 



20 



K = -^hr F { c ^) d ^ + c{ V )drf + c( v f )drf}, (52) 

Zi COo (7 J/J/ 

where the C(P) functions are given by, 

C(rj) = (sin0 + V2cos0) , 

C(rf) = (V2 sm9 -cos fl) . (53) 

The hadronisation into one pseudoscalar meson P of the quark bilinears in T^o proceed via 
the pseudoscalar currents P l . Concretely, P°, P 3 and P 8 , whose expressions in terms of one 
P meson can be obtained from Eq. (26). This leads to: 

P 3 = 2B Ftt°, 

P 8 = 2B F (cos Ori + smOri') , 

P° = 2P O -F(-sin0?7 + cos0?/) . (54) 
Finally, by putting all together, we get the following result for the branching ratio: 

BR(r - MP) = !■ AV2 (<<^) I V | T | 2 , (55) 

hf 

where the X(x, y,z) function is defined in Eq. (43) and again T T is the total decay width of 
the r lepton. The averaged squared amplitude is given by, 

\ J2 l T ' 2 = T~ J2 [ 2m » m r (4°p* - VpVp*) + Or +ml~ m 2 P ) (a k P a™* + b k P b™*)] , 

i,f k,m 

(56) 

with k,m = Z, A , and 

4 = - 9 — a ^^P(m r -m M )(P L + P fi ) , 
2 cosfe'iy 2 mi 

1,2 



2 cos t% 2 



a 9 

(In = 



- - ^ ^ 0, (") K' + <*'). 

^ - ^4(^>-<^(^-*n- (57) 

The B^k(P) functions are given, correspondingly, by the third entry of: 



,2 / „(p)* „(p)* 

B^(n) = - 
'■ 11 sin/3 cos/3 
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B?(v) 



4^3 
1 



-a. 



(p)* 



sin (3 



—a, 



(p)* 



sin /3 



2 2 



2 (sin# + v^cos^) +^-^ (3m^-4m^)sin^ + 2v / 2m^cos^ 

/ COS O L 



(p)* 

(P)* r 



/ \ o- w * r i 

( cos6> - y/2sme) + ^—5 (3m 2 - Am 2 K ) cos9 - 2\f2m 2 K sin6> 

\ / COS £5 L 



Bt\P) = B^*(P), (58) 

where the af\ functions are defined in Eq. (36). Notice that in this Eq. (58) the relations 
between the quark and meson masses of Eq. (18) have been used again. 



3.3 Predictions for r — > jip and r — ► 



The r — > pp° decay is related to the r — > pn + 7i^ channel since the p decay proceeds mainly 
to 7r + 7r~. Indeed a p° is not an asymptotic state : the experiment reconstructs its structure 
from the two observed pions. In addition, from the chiral point of view, two pions in a 
J = I = 1 state are indistinguishable from a p. Therefore one has to define the branching 
ratio of r — > pip in close relation to that of r — > pn + -K~ as follows: 

^max 

BR(r^/ip°) 



1 



64 7r 3 m 2 . T r 



ax Amax 
min *J ^min ^ J* 



(59) 



where T 7 is defined in Eq. (37) and all functions and form factors involved are as those of 
t — > pn + 7i~ decay, with the exception of the integration limits in s which are now: 



M 2 p - \m p Y p , 



1 



W = M p 2 + -M p r p . (60) 

Similarly, the r — > pip decay is related to the r — > pK + K~ and r — > pK°K° decays since 
the decays proceeds mainly to K + K~ and to K°K°. Therefore, we define: 



BR(r -> p0) = 



64 7r 3 m 2 T r 




1 



jst+.r:- 



+ 



i,/ 



(61) 



where again T 7 is defined in Eq. (37) and all functions and form factors involved are as 
those of r — > pK + K~ and r — > pK°K° correspondingly, except for the integration limits in 
s which are now: 



Ml - -MJ\ 



Ml + -M,r d 



(62) 
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In Eqs. (60,62), T p = T p (M 2 p ) and I> = T (M|). 

4 Numerical results and discussion 

In this section we present the numerical results of the LFV semileptonic r — > fiPP and 
t — > \xP decay rates within the constrained MSSM-seesaw scenarios described in Section 2. 
Since our main goal is to explore if the predicted rates can or cannot reach the present 
experimental sensitivities we will focus mainly on choices of the input parameter values that 
lead to large £32 and therefore to large LFV semileptonic r decay rates. As we have seen 
in the previous Section 2, within the scenario with hierarchical heavy neutrinos and for 
$1,3 = 0) this means large values of 82 and large values of tojv 3 - On the other hand, since 
all these rates grow with tan /3, in the following numerical analysis we will focus mainly on 
large tan (3 values. In the first subsection we will present the numerical results, from our full 
computation of the LFV semileptonic tau decay rates and will explore the dependence with 
the most relevant parameters in the constrained MSSM scenarios. In the second subsection 
we will include a comparison between our full and some approximate results in the large 
tan p region. Moreover, we will also analyse to what extent the Higgs dominance hypothesis 
holds for these LFV semileptonic r decays and compare our predictions with other results 
in the literature. We will conclude this section by showing that for some particular choices 
of the input parameters, the rates for some channels indeed reach the present experimental 
sensitivity. 

4.1 LFV semileptonic tau decay rates 

Firstly, we present the results for the simplest case of 9 2 = and study the relative im- 
portance of the various contributions to the decay rates that have been presented in the 
previous section. Then we explore the increase in the rates for larger values of 9 2 . Since we 
are setting in the whole numerical analysis A = and sign (//) = +1, the relevant SUSY 
parameter will be M S usy = M = Mi/ 2 . In the study of the behaviour of the rates with 
Msusy we pay special attention to the decoupling or non-decoupling behaviour of the SUSY 
particles at large M SUSY . 

In Fig. 4 we display the prediction of BR(r -> fiPP), with PP = tt+tt-, K+K~, K°K°, 
7r°7r°, as a function of M SU sy and for the particular choice of 0j = (i — 1, 2, 3). We consider 
both CMSSM (left panel) and NUHM (right panel) scenarios. We set here tan (3 = 50 and 
our "reference" values of m Nl 2 :i = (10 10 , 10 11 , 10 14 ) GeV. For the NUHM case we set in 
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Figure 4: BR(r -> /j,PP) for PP = Tr+rr" , K+K~ , K°K° , 7r°7r° as a function of M SU sy = M = M 1/2 in the 
contrained MSSM-seesaw scenarios: CMSSM (left panel) and NUHM (right panel). 



addition 8\ = —1.8 and 62 = 0, which have been shown in [24] to lead to low Higgs boson 
mass values. Concretely, for Q { = and 250 GeV < M SU sy < 900 GeV the predicted masses 
are within the range 110 GeV < m A o,m H o < 180 GeV, which are indeed very close to their 
present experimental bounds. 

The first obvious conclusion from Fig. 4 is that the rates of the different channels exhibit 
the following hierarchy, BR(r -> ^tt+tt-) > BR(r -> jj,K + K~) > BR(r -> fiK°K°) > 
BR(r — > fi7T°7T ). This hierarchy can be understood in terms of the dominant electromagnetic 
contribution and the relative phase space suppression. We also see that the decoupling 
behaviour for large Msusy is clearly manifest in the universal case, where all the rates 
decrease as Msusy grows. In contrast, it turns out that, in the NUHM case, the decoupling 
behaviour is only manifest in the r — > /i7r + 7r~ channel, whereas the r — > fiK + K~, r — > 
fiK°K° and r — > fiTT°n rates do not decrease with Msusy in the large Msusy region. This 
behaviour can be better comprehended by analysing separately the different contributions 
to these channels, as shown in Fig. 5. 

The results displayed in Fig 5 for the r — ► yU7r + 7r~ channel show the dominance of the 
photon-mediated contribution in this case, which is in fact indistinguishable from the total 
rate in this plot, for all the explored parameter values. The Higgs- mediated contribution is 
subdominant by far due to the highly suppressed couplings of the Higgs to the light u and d 
quarks, which after the hadronisation of the corresponding quark bilinears result in Htt + tt~ 
couplings proportional to (see Eq. (41)). This plot also exhibits the non-decoupling 
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Figure 5: Rates of the various contributions to BR(t — > fiK + K ) (upper left panel), BR(r — > /j,K°K°) 
(upper right panel), BR(r — > fiir + Tr~) (lower left panel) and BR(t — * /irj) (lower right panel) as a function 
of Msusy = Mo = M1/2 in the NUHM scenario. 



behaviour of the SUSY particles in the Higgs-mediated contribution. The particular pattern 
of this contribution as a function of Msusy is a consequence of two facts. First, the well 
known constant behaviour with MgusY of the LFV Hr[i form factor at large Msusy- Second, 
the encountered Higgs mass behaviour with Msusy ; analysed in [24], which, for this choice 
of 8^2 and for the studied Msusy interval, first grows softly, reaches a maximum and then 
decreases softly. 

The r — > / U7r°7r° channel is only mediated by the Higgs bosons and a similar suppression 
of Htt tt couplings as in the H-k + tt~ case occurs, leading to very low predicted rates. These 
low rates and the non-decoupling behaviour of this channel can be clearly seen in Fig. 4. 
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Figure 6: BR(r — > /iP) for P = zr, 77, 77', p, as a function of Msusy 
MSSM-seesaw scenarios: CMSSM (left panel) and NUHM (right panel). 



Mq — M1/2 i n the contrained 



The results for the r — > fiK + K~ channel that are depicted in Fig 5 are interesting 
because the photon- and the Higgs-mediated contributions compete in this decay. In fact 
the Higgs-mediated contribution can equalise, or even exceed that of the photon, dominating 
the total rate in the large Msusy region. Both photon- and Higgs-mediated contributions 
are similar around Msusy = 750 GeV. The reason for this larger Higgs contributions than in 
the previously studied tttt case is because of the larger Higgs couplings to the strange quarks 
which result in HKK couplings proportional to m 2 K (see Eq. (41)). 

The results for the r — > fiK°K° channel in Fig 5 are very similar to those for r — > 
/j,K + K~. One difference is the point where the Higgs-mediated contribution crosses the 
photon one, which for r — > fiK°K° is around Msusy = 700 GeV. Another interesting 
difference is that this rate is always slightly smaller than r — > fiK + K~ due to the fact 
that the photon-mediated contribution to r — > fiK°K° occurs just by the meson resonances, 
whereas the r — > fxK + K~ channel can also be mediated via pure electromagnetic interaction. 
This difference is clearly summarised in the several contributions to the Fy +K and Fy°^° 
form factors in Eq. (B.3) of Appendix B. 

The predictions of BR(r — > ^P), with P being here a pseudoscalar meson 7r,T],T]' or 
a vector resonance p, 0, as a function of M S usy are displayed in Fig. 6. We also consider 
CMSSM (left panel) and NUHM (right panel) scenarios. In the universal case we find 
the following hierarchy, BR(r — > pp) > BR(r — > fief)) > BR(r — > pr]') > BR(r — > prj) 
> BR(r — > /i7r). We obtain again the expected decoupling behaviour for large M S usy m 
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Figure 7: Light m^o predictions as a function of non-universal parameters Si and 82 in the NUHM scenario. 
The predictions for m#o are indistinguisable from those of m^o in this figure. 



this universal scenario, while in the NUHM scenario the non-decoupling behaviour is clearly 
manifest for r — > fit], r — > firj' and t —>■ fin. The r — > /ip rates in the NUHM scenario are 
the largest ones, except in the large Msusy region, where r — > \n] and r — > fir]' rates exceed 
them. These two channels are by far dominated by the Higgs-mediated contributions in the 
full MgusY explored interval, as can be seen for the r] case in Fig 5. The reason for this 
Higgs dominance is because of the large Higgs couplings to the strange components of the 
r) and 7/ mesons, which result in large A — rj and A — rj' "mixings" proportional to m 2 K as 
explicitely given in Eq. (58). 

One of the most important outcomes from the previous analysis, corresponding to the 
9i = choice, is that for both scenarios and for the chosen input parameters, the predicted 
rates for both r — > fiPP and r — ► fiP channels do no reach their corresponding experimental 
bounds, and even in the best cases of r — > fi7f + 7t~ and r — * fip they are still two orders of 
magnitude below their present experimental sensitivities. In the following, we will therefore 
focus on larger values of 9 2 - 

In order to reach the larger rates as possible in the 0$ 7^ case, one needs to explore 
first the optimal values of S\ and 62 which lead to light Higgs bosons. We summarise the 
predictions for the relevant Higgs boson mass, m^o (and m#o), as a function of Si and S2 
in Fig. 7 for the extreme value of 82 = 2.9e J7r / 4 . The reason for this particular choice of 
62 is due to the fact that it leads to the maximum value of \S^2\ which is compatible with 
our hypothesis of perturbativity, as shown in Fig. 1 and discussed in Section 2. We have 
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Figure 8: Light m^a predictions as a function of Msusy = Mq = M 1 / 2 in the NUHM scenario. The 
predictions for m#o are indistinguisable from those of m^o in this figure. The predictions for the CMSSM 
scenario (Si = S2 = 0) are also included for comparison. 



chosen here tan/3 = 50 and two representative values of Msusy = 250 and 750 GeV for 
moderate and heavy SUSY spectra, respectively. The other parameters are set to the values 
of m Ni = (10 10 , 10 11 , 10 14 ) GeV, Q x = 6 3 = 0, A = and sign(» = +1. To ensure that our 
results are indeed experimentally viable, we have included in this, and in the following figures, 
only the solutions where the three neutral Higgs boson masses are above the experimental 
bound for the lightest MSSM Higgs boson, which at present is 110 GeV for tan/3 > 5 (99.7% 
C.L.) [7]. The most interesting solutions with important phenomenological implications are 
found for negative 5\ within the range (—3, —2) and very small and positive 82, the choices 
selected for Fig. 7. In this figure, for all the explored values of Si and 62 , we find a value of 
m^o that is significantly smaller than what one would encounter in the universal case (here 
represented by the choice Si = 62 = 0). This is truly remarkable in the case of large soft 
breaking masses, as can be seen, for instance, in the panel with Msusy = 750 GeV, where 
low values of m^o ~ 150 GeV are still found. 

The behaviour of the predicted m A o as a function of Msusy is depicted in Fig. 8. Here the 
specific values of Si = {—2.45, —2.4, —2.35, —2, —1, 0} and S2 = 0, 0.2 have been considered. 
This figure illustrates again the interesting departure in NUHM scenarios from the linear 
behaviour of m^o with Msusy, which is generic in the universal case (Si = S2 = 0). The same 
pattern with M S usy was also found for the 9^ = case in [24], but obviously for different 
choices of Si and 82- 

The corresponding predictions for 9 2 = 2.9e l7f//4 of the nine LFV semileptonic r decays 
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Figure 9: Predictions of BR(r — > I^PP) and BR(r — > /xP) as a function of Msusy in the NUHM scenario 
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given in Table 1. 



studied in this work as a function of Msusy are shown in Fig. 9. In this case, we work 
with 5i = —2.4 and 62 = 0.2, that drive us to Higgs boson masses around 150 GeV even 
for heavy SUSY spectra, as can be seen in Figs. 7 and 8. In this Fig. 9 we can see that, 
compared to predictions in Figs. 4 and 6, the new choice of 82 increase all the rates about 
two orders of magnitude. All the rates exhibit the same hierarchy as in the previous plots, 
being BR(r — ► /i7r + 7r - ) and BR(r — ► up) the largest ones. Indeed, the predictions of these 
two latter channels reach their present experimental sensitivities at the low Msusy region, 
below 200 GeV and 250 GeV respectively, for this particular choice of input parameters. 

4.2 Comparison between the full and approximate results 

It is interesting and useful to provide simple formulas that can approximate reasonably well 
our full predictions. The most popular approximation when predicting LFV rates is to work 
with expressions that are valid only in the large tan/? region. The justification for this is 
obvious since all these LFV rates are known to grow with tan {3. It is specially important in 
scenarios where the LFV rates are dominated by the Higgs mediated diagrams, since these 
latter grow much faster with tan /3 than the photon or Z boson mediated ones. Accordingly, 
we will pay more attention to the semileptonic r — > pP and r — ► pPP channels that can 
be dominated by the Higgs bosons and whose present experimental sensitivities are the 
best ones. This leads us mainly to the r —>■ prj and r — > pK + K~ channels. The other 
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approximation which is used frequently in the literature, due to its simplicity, is the mass 
insertion approximation, where the tau-muon LFV is encoded in the dimesionless parameter 
#32, already introduced in Sec. 2, and whose expression in the LLog approximation is given 
in Eqs. (9) and (10). 

We start by considering the large tan (3 limit of the tau-muon-Higgs form factors that 
are the relevant ones for the LFV Higgs-mediated processes. The full one-loop Higgs form 
factors were computed in [52] (see also [21]) and are collected in Appendix A. 2. At large 
tan/3, H L dominates H R by about a factor m T /m^. Moreover, and Hf° are by far the 
largest form factors in this limit, and one can safely neglect . More specifically, by using 
the mass insertion approximation, their chargino and neutralino contributions in the large 
tan/3 limit give, correspondingly, the following results [52], 

< = l < = ^l2^: fetan/3 ' (63) 

<> = ,< 0) = ^^(l-3tan 2 M 53 2 tan^. (64) 

One can further verify that H c dominates H n by about a factor 20, so that in the following 
we will take Hl — Hl, c - 

On the other hand, we also consider the large tan f3 limit of the functions that define the 
Higgs couplings to one meson, B(P) in Eq. (58), and to two mesons, J(PP) in Eq. (41). It 
leads to the following results: 

i-^=tan/3 (3ml — Am 2 K ) cos 9 — 2\p2m\ sin#J , 
(3m\ - Am 2 K ) sin 9 + 2\/2m 2 K cos o] , 
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By using the above sequence of approximations and by neglecting the muon mass, we finally 
get the following simple results, 
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where s max , s min , t max and t min are given in Eq. (43). The results for the other channels 
can be similarly obtained by using the corresponding B(P) or J(PP) functions and the 
corresponding meson masses (with an additional 1/2 factor in the case of BR(t — > fiir°iT ) 
to account for identical final state particles). We get, 

BR(, = 1.5 x 10- fel 2 (^^y (^) 6 . (68) 
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In all the above approximate results of the LFV semileptonic tau decay rates we see ex- 
plicitely the strong dependence with both tan/5 and the corresponding Higgs boson mass, 
being (tan/3) 6 and (1/m^) 4 , respectively, which are characteristic of the Higgs mediated 
processes. 

Regarding the comparison with other works, first, we notice that our numerical prediction 
for BR(r — > /irj) in Eq. (66) does not agree with the original estimate in [31] that gives a 
decay rate a factor 7 larger than ours. We believe that the discrepancy comes from our 
different approaches to describe the hadronisation of quark bilinears. Our numerical result 
is closer to that in [32] whose prediction is larger than ours in a factor of 2. Notice, that the 
comparison with this latter work must be done by switching off the bottom-loop induced 
contributions and the higher order loop-effects enhanced by tan (5 factors which were taken 
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Figure 10: Comparison between the predicted rates of BR(r — > firf) in the NUHM from our full 1-loop 
computation and from the approximate result of Eq. (66) as a function of tan (3 (left panel) and AfsusY = 
Mo = Mi/2 (right panel). The horizontal lines are the present experimental bounds given in Table 1. 



into account in [32] but we are not including here. This means setting their £ q parameters 
to £b = and £ s = 1 in their formulas. We believe that this small discrepancy is mainly due 
to the different approaches for hadronisation. In particular, they neglect the m u ^ masses 
whereas we are taking into account chiral symmetry breaking effects via the explicit and 
m 2 K dependences, which are well determined in the ;^PT approach. On the other hand, our 
prediction for BR(r — > fir]') is slightly above BR(r — > firj), due basically to the larger Higgs 
coupling to r/, \B L (rf)^\ > \B L {r])^ Ao) \. The prediction in [32] of BR(r — > fir)') is, however, 
a factor 100 smaller than ours. The prediction for BR(r — > fin) here and in [32] agree within 
a factor of 2. Finally, the prediction for BR(r — ► fiK + K~) in [34] is larger than our result 
in about a factor 50. 

The goodness of the above approximate result for r — >• firj in Eq. (66) can be seen in 
Fig. 10, where it is compared with the full result as a function of tan/? and Msusy- It is clear 
that, for tan/5 values larger than about 30 the approximation is quite good, providing rates 
that are at most a factor of 2 above the full predictions. Moreover, the behaviour with tan (3 
of the full result at this region is well described by the (tan f3) 6 behaviour of the approximate 
one. Regarding the behaviour with Msusy, we see again that the approximate and full 
results differ by less than a factor of 2 and they both follow the same pattern. The displayed 
dependence with Msusy can be easily understood from the dependence of m^o with this 
parameter, as was shown in Fig. 8. For the studied range in this plot, 250 < Msusy (GeV) < 
650, this leads to a relatively small variation in the rates of about BR max /BR min ~ 5. 
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The Higgs dominance approach, however, is not so good for other LFV tau decay chan- 
nels. In particular, it is clearly not a good approximation for t —>■ 3fi because, in this case, 
there are other contributions from 7-mediated, Z-mediated and box diagrams that enter 
into the full computation [2,21]. In the NUHM scenarios that are considered here with 
small Higgs masses, one may guess that the Higgs mediated contribution could dominate 
the rates at large tan/3, but it is not so as will be shown next. By performing a similar 
analysis as we have done before, that is, by using the tau-muon- Higgs form factors in Eq. 63 
and plugging it into the exact formula for the Higgs-contribution [21], we get in the large 
tan/3 limit, 
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which is in good agreement with the original result in [18] and also with posterior esti- 
mates [19,32]. The comparison between the full (i.e. including one- loop SUSY diagrams 
mediated by 7, Z, ho, Hq, Aq and box diagrams which are taken from [21]) and the ap- 
proximate numerical results for this channel is shown in Fig. 11. We see that the formula 
in Eq. (73) predicts rates that are about a factor of 2 larger than the exact Higgs-mediated 
contribution. Therefore, for large tan/3 values, it provides a good estimate of the Higgs 
contribution. However, the total rates are much larger than the Higgs contribution, since 
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the photon-mediated diagrams give by far the dominant contribution in this channel. For 
instance, we see in Fig. 11 that the total and Higgs rates differ in about two orders of mag- 
nitude for Msusy ~ 250 GeV and in more than a factor 5 for M S usy ~ 750 GeV. It is 
remarkable that, in this channel, the photon dominance holds largely even in scenarios with 
very heavy SUSY spectra, as for M guSY ~ 750 GeV, and Higgs bosons as light as m H = 160 
GeV. Therefore, the total rates for this channel can be better approximated by the simplified 
formula of the photon-mediated contribution, 

BR(r^3,) 7 ._ = iL( log ^-ii)BR(^,n) 

= 2.3 x 10~ 3 BR(r -> // 7 ) 

o , 1ft .5, r ,2/ 100 V /tan/A 2 

where the last line has been obtained by using the result of BR(r — > ^7) in the mass insertion 
approximation for equal SUSY mass scales and in the large tan/3 limit. It is also interesting 
to compare this estimate with the present experimental upper bound for this channel which 
is 3.2 x 10~ 8 [16,17]. We see in Fig. 11 that, for the chosen parameters in this plot, the 
predicted rates are above the present experimental bound for Msusy < 300 GeV. 

Similarly to the r — > 3/x channel, the semileptonic r — > fiPP decays (with the exception 
of t — > /i7T 7r°) are clearly dominated by the photon contribution and therefore they can 
be better approximated by the corresponding simplified formulas of this contribution. By 
neglecting the \i mass we have found the following approximate result, 



x ihf^^^ (76 > 



And from this we get, 

BR(r - /x7r + 7r-) 7approx = 2.5 x 10" 3 BR(r — > ^7) 



100 \7 tan ^ 2 



BR(r — ► fJ-K + K~) lapprox = 2.0 x 10~ 4 BR(r - m ) 

v MsusY(GeV) 
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Figure 12: Comparison between the full rates and the approximate results in the NUHM scenario by con- 
sidering just the photon-mediated contributions for r — > /i7r + 7r~ (left panel) and r — > fiK + K~ (right panel) 
as a function of Msusy = Mq = Mi/ 2 . The dashed horizontal lines are the present experimental upper 
bounds [25]. 



BR(r - »K°K% wiOK = 1.2 x 1(T 4 BR(r ^ mt) 



BR(r — > /ip) 7approx = 2.3 x 1(T 3 BR(r - /i 7 ) 



BR(r - /i0) 7approx = 8.4 x 10~ 5 BR(r — > fij) 

stc .2 / 100 \ 4 /tan /A 2 , . 

As can be clearly seen in Fig. 12 these results approach pretty well the full rates for most 
of the M SU sY studied region. For BR(r — > fi7r + ir~), they are indeed indistinguishable in this 
plot. It is only at very large Msusy > 750 GeV that the approximate result of BR(r — ► 
/j,K + K~) separates slightly from the full result, due to the Higgs contribution which competes 
with the photon one in this region. 

For completeness and comparison, we also include here the predictions for the leading 
LFV tau decay channel, r — > /ry. Fig. 13 displays the predictions of the full and approximate 
rates for this r — > p>q channel. The full rates are taken from [21] and the approximate ones 
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Figure 13: Comparison between the full 1-loop prediction [2] and approximate results of Eq. ( 82) for r — > fij 
as a function of Msusy = Mq = Afi/ 2 in the NUHM scenario. The horizontal line is the present experimental 
upper bound [15]. 



are given by the result of the MI approach [3,33], which at large tan/3 and for equal SUSY 
mass scales is, 

BR(r — > /i7) approx = ttttz 2r . Aa ., 4 tan p) 

144007H T T sin 9 W Mg USY 

In this case, and for the chosen parameters in this plot, the approximate and the full results 
agree to better than a factor 2. We have verified, however, that for other choices of Si^ the 
difference between them can be larger. Regarding this difference, we emphasise that in using 
the MI approach and LLog approximation one has to be carefull because they are known to 
fail in some regions of the CMSSM parameter space. For instance, in [53], the departure of 
the MI from the exact result is estimated to be up to 50% for |5 32 | ~ 1. In [22] it has been 
found that the use of the MI and LLog for large trilinear couplings, Aq ~ O (1 TeV), can 
fail in several orders of magnitude. 

The most evident conclusion from Fig. 13 is that for the chosen parameters in this plot 
and for Msusy < 1600 GeV, the r — > /ry rates are above the present experimental sensitivity, 
therefore this tau decay channel is at present the most competitive one in setting bounds on 
the tau-muon LFV. However, besides experimental issues, the limitation of this channel is 
that it is not sensitive at all to the Higgs sector. In this sense, the semileptonic channels are 
more interesting, and can be clearly competitive in the large Msusy ~ 0(1 — 2 TeV) region. 
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In Figs. 14 and 15 we plot finally the predictions for BR(r — > ^,K + K~) and BR(r — > /j,rj) 
as a function of one the most relevant parameters for these Higgs-mediated processes which 
is the corresponding Higgs boson mass. Firstly, we see again that the approximate and 
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Figure 14: Predictions for BR(t — > /iK + K~) as a function of m#o in the NUHM scenario. A comparison 
between the full 1-loop computation and the approximation given by Eq. ( 67) for various choices of large 
tan/3 and mj\r 3 is included. The horizontal line is the present experimental upper bound [25]. 



exact results of the Higgs contribution agree within a factor of two for both channels, but 
the agreement of the full result with respect to the Higgs contribution is clearly worse in the 
case of t — > fiK + K~ than in r — > \ir\. In the latter, the agreement is quite good because 
the Z-mediated contribution is negligible, and this holds for all Msusy values in the studied 
interval, 250 GeV < M SU sy < 750 GeV . In the first, it is only for large M SUS y that the 
if-mediated contribution competes with the 7-mediated one and the Higgs rates approach 
the total rates. For instance, Fig. 14 shows that for M SU sy = 750 GeV and m#o = 160 GeV 
the total rate is about a factor 2 above the Higgs rate, but for m#o = 240 GeV it is already 
more than a factor 5 above. 

In these figures we have also explored larger values of rriN 3 and tan/3, by using in those 
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Figure 15: Predictions for BR(r — > /irj) as a function of m^o in the NUHM scenario. A comparison between 
the full 1-loop computation and the approximation given by Eq. ( 67) for various choices of large tan /3 and 
mjv 3 is included. The horizontal line is the present experimental upper bound [15]. 



cases the approximate formula, and in order to conclude about the values that predict 
rates comparable with the present experimental sensitivity. We can conclude then that, at 
present, it is certainly r — > fit] the most competitive LFV semileptonic tau decay channel. 
The paremeter values that provide rates being comparable to the present sensitivities in this 
channel are tan f3 = 60 and mjv 3 = 10 15 GeV which correspond to [5^1 — 2. These large rates, 
however, should be taken with care and be considered just as an order of magnitude estimate 
since, as we have explained in Sec. 2.1, they correspond to neutrino Yukawa couplings which 
are clearly in the non-perturbative regime. This is why we do not provide the corresponding 
full rates for them. 
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5 Conclusions 



In this paper we have presented a complete one-loop computation of the branching ratios 
for the LFV semileptonic r decays within the context of two constrained MSSM-seesaw 
scenarios, the CMSSM and the NUHM. We have included both analytical and numerical 
results for the particular channels: r — > fiPP, with PP = 7r + 7r~, 7r°7r°, K + K~ , K°K°; r — > 
fiP with P = n,r], r)'\ and r — > /xp, r — > The analysis of the channels r — > fiPP, 

with PP = 7r + 7r~, 7r°7r°, K°K°, and r — > //p, r — > p0 are, to our knowledge, the first ones 
in the literature within the CMSSM-seesaw context. In addition, we have compared our 
predictions for r — > fiK + K~ and for r — > pP with P — 7T,r), f]' with previous predictions in 
the literature and found some discrepancies. 

Our treatment of hadronisation has involved two different procedures : for the 7 con- 
tribution we have employed state of the art form factors (detailed in Appendix B), as this 
amplitude is fairly dominated by resonance states; for heavier intermediate contributions (Z 
and Higgses) we have a local (point-like) vertex driven by chiral symmetry. It is difficult 
to estimate the errors of this procedure. For the hadronization of the 7, large- Nq inspired, 
the error should be smaller than 30% [46] (at amplitude level), based on the fact that sub- 
leading terms in the expansion have been included through the widths of resonances. The 
hadronization of currents driven by the Z or Higgses, on the other side, is less known and it 
is not possible to give a reliable error estimate. 

Our results for r — > p7r + 7r~ demonstrate that this channel is clearly dominated by the 
photon-mediated contribution in all the studied region of 100 GeV < M SUSY < 1000 GeV. 
In fact it is by far, the r — ■> pPP channel with the highest rates, reaching values close to its 
present experimental bound at 4.8 x 10 -7 for some input parameter values. Concretely, it 
happens for low M SUSY ~ 100-200 GeV, large tan/3 ~ 50-60, large m Na ~ 10 14 -10 15 GeV 
and large arg(# 2 ) ~7r/4 — tt/2 (these two latter parameters producing a large £32 ~ 0(1)). 
In contrast, r — > fnr°7r° can only be mediated by h° and H° Higgs bosons and their rates 
are very small. Besides, they are not yet comparable with data, since there is no bound 
in this channel. The cases of r — > fiK + K~ and r — > fiK°K° decays, are much more 
interesting. In these two channels, the photon-mediated contribution dominates in most of 
the studied region of Msusy, except at large, Mgugy > 750 GeV values, where the Higgs- 
mediated and the 7-mediated contributions can compete. This competition happens in 
specific constrained scenarios of NUHM type with low m#o ~ 100 — 200 GeV values and 
very heavy SUSY spectrum with MgusY > 750 GeV. This peculiar MSSM spectrum and 
the fact that Higgs bosons couple stronger to K + K~ (and K°K°) than to 7y + ty~ (and 7r°7r°) 
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is the reason why the H- and 7-mediated contributions can compete in r — > fiK + K~ but 
not in r — > fnv + n^ . Furthermore, due to the fact that the photon diagram still dominates 
BR(r — > fiK + K~) in a large region of the parameter space with 100 GeV < M S usy < 750 
GeV, the involved hadronic form factors do play a crutial role in the final rates. Consequently 
our results for this channel are in disagreement with those of [34] where they only included 
the Higgs-mediated contribution. We have also shown that the largest predicted rates for 
BR(r — > fiK + K~) are, as in r — > fj,ir + ir~, at the region with low M S usy ~ 100 — 200 GeV, 
large tan (3 ~ 50 — 60, large m N3 ~ 10 14 — 10 15 GeV and large arg(# 2 ) ~ tt/4 — ir/2 values. 
However, the predicted rates do not reach yet the present experimental sensitivity, which in 
this channel is at 8 x 10~ 7 . 

Our results for r — > \ir) and r — > fir/' demonstrate that these two channels are largely 
dominated by the A°-mediated contribution and their predicted rates are very competitive in 
the case of NUHM scenarios with low m A o ~ 100 — 200 GeV values and large tan (3 ~ 50 — 60. 
This is in qualitative agreement with previous estimates in the literature. However, we 
have found some important numerical discrepancies with respect to the estimate in [31]. 
Concretely, the predicted rates in the present work are smaller than those in [31] by a 
factor of about 7. We believe that these discrepancies are due to the different procedures 
of quark bilinear hadronisation. We claim that our results which are based on the well 
defined and more refined hadronisation prescription by provide a better estimate. The 
rates for BR(r — > firj) have also been compared with those in [32,33] which are within 
the different context of non-constrained MSSM and with input £32 not being connected to 
neutrino physics nor seesaw mechanism. We have checked, that the predicted rates are in 
reasonable agreement with these two works for, £32 ~ 0(1) , which in our case is reached by 
input seesaw parameters of mjv 3 ~ 10 14 — 10 15 GeV and large arg(# 2 ) ~ 7r/4 — 7r/2. 

In addition, we have presented in this work a set of useful approximate formulas for all the 
semileptonic r decays that we have compared with the full-one loop results and concluded 
that they give reasonable good estimates, say differing in less than a factor of two respect 
to the full result. We have also compared these results with those for the leptonic channel, 
r — > 3/x, and the radiative decay, r — > /ry. 

Our overall conclusion is that, for the same Constrained MSSM-Seesaw input parameters, 
r — >• /ry is the most competitive r decay channel in testing the values of the LFV parameter 
£32, but it is not sensitive at all to the Higgs sector. Interestingly, the most competitive 
channels to explore simultaneously LFV and the Higgs sector are r — > fir], r — > /irj' and also 
t — > fiK + K~ . The t — > fiK + K~ channel is certainly more efficient than r — > 3 \x as far as 
the sensitivity to the Higgs sector is concerned. Otherwise, the golden channels to tackle 
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the Higgs sector are undoubtly r — > /irj and r — > //?/. On the other hand, the rest of the 
studied semileptonic channels, r — > fi7v + 7i~ , etc., will not provide additional information on 
LFV with respect to that provided by r — > /ry. 
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A LFV form factors 

In this Appendix we collect the main analytical formulae containing the full 1-loop results 
of the SUSY contributions to the relevant r — fx LFV form factors for the present work, cor- 
responding to: 7r/i, Zrji and Ht\i vertices. All the couplings and loop functions appearing 
in the following formulae are defined in [21,24]. 

A.l Form factors for the 777/ vertex 

Our convention for the form factors A\^ defining the 777/. vertex is as follows: 



where q is the off-shell photon momentum, P LjR = (1 =f 7s)/2, e is the electromagnetic 
positron charge and m T is the r lepton mass. 

In the SUSY-seesaw context there are one-loop contributions to these form factors that 
come from the chargino and neutralino sectors respectively, 



te [q 2 la (A^P L + A*P R ) + im T <j al3 q p {A L a P L + A*P R )] , 



(A.l) 



L,R _ j^(n)L.R _|_ j^[c)L,R 



a = 1,2. 



(A.2) 



The neutralino contributions are given by, 
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Figure A.l: Relevant SUSY one-loop diagrams for the photon-mediated contributions to LFV semileptonic 
r decays. 
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where xax = /mf and the indices are A = 1, .., 4, X = 1, .., 6. 
The chargino contributions are given by 
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A(c)R = a(c)L\ 



(A.8) 



where in this case xax = m ^-/ m o x an d the indices are A = 1,2, X = 1,2,3. Notice that 
in both neutralino and chargino contributions a summation over the indices A and X is 
understood. 



A. 2 Form factors for the Zt\x vertex 

Our convention for the form factors F L R defining the Zt\i vertex is as follows: 

[FlPl + F R P R ] . 



(A.9) 
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Figure A. 2: Relevant SUSY one-loop diagrams for the Z-mediatcd contributions to LFV semileptonic r 
decays. 



The Z-boson form factors have also the two kinds of contributions, from neutralinos (n) 
and charginos (c), 

Fl {R) = F^ R) + F$ R) . (A.10) 
The results for the corresponding form factors are the following: 
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where again the indices are A,B — 1, ..,4, X, Y = 1,..,6 in the contributions from the 
neutralino sector and A, B = 1,2, X, Y = 1,2,3 in the contributions from the chargino 
sector, and a summation over the various indices is understood. 



A. 3 Form factors for the Hrfi vertex 



Our convention for the form factors H^ R defining the H p Tfi vertex is as follows: 
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Figure A. 3: Relevant SUSY one-loop diagrams for the Higgs-mediated contributions to LFV semileptonic r 
decays. 



As in the previous cases, we separate the contributions from the neutralino and chargino 

sectors, 
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The results for the form factors are the following, 
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Correspondingly, the result for the chargino contribution H^iR) c can De obtained from the 
previous H^Ir) n D Y replacing everywhere, 
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In the previous formulae, the index p refers to the each of the Higgs bosons. Concretely, 
Hp = h°,H°,A° for p = 1,2,3, respectively. The other indices are again A, B = 1,..,4, 
X, Y = 1, .., 6 in the contributions from the neutralino sector and A, B = 1, 2 and X, Y — 
1, 2, 3 in the contributions from the chargino sector. A summation over all the indices is also 
understood. 
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B Hadronic form factors 



Our construction of the vector form factors Fy P (s), defined by Eq. (28), follows the idea 
put forward in [54] that lie on two key points : 

1/ At s <C M\ (being Mr a generic resonance mass), the vector form factor should match 
the 0{p A ) result of %PT. Hence our form factors will satisfy the chiral constraint. 

2/ Form factors of QCD currents should behave softly at high transfer of momenta [38], 
i.e. they should vanish for s ^> M\. Accordingly we will demand to our form factors 
that they satisfy this asymptotic constraint. 



In the Nq — > oo limit resonances have zero- width. However those present in the relevant form 
factors in tau decays do indeed resonate due to the available phase space. As a consequence 
we need to include energy-dependent widths for the wider resonances p(770) and p(1450), 
or constant for the narrow ones : cj(782) and 0(1020). For the p(770) we take the definition 
put forward in [55] : 
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The 0(p A ) determination of the vector form factors was done in [49]. Requiring that our 
expressions match that result at small transfer of momentum we get the following expres- 
sions : 
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8111 ° V M*-3-iM„T u + C ° S 9V Ml-s-iM,Y 4 



x exp 3 Re ^Hkk(s) j 



7 K°K° 
V 



M 2 



2 M 2 - s - iM p V p {s) 



exp 



+ 



sin 2 # 



Ml 



v 



Ml-s- iM u T u 



2Re(H mr (sj) + Re(H KK (s)) 

- + COS Uy 



Mi-S- iM.Xr 



x exp 



3 Re [Hkk{s) 



where we have used the definitions : 

0, 



P = 



7 
6 

Hpp(s) 
M P (s) 

Ms) 

k P ((j,) 



3M 2 ' 

FyGy 

F 2 

FyGy 



F 2 



(1+/3) - 1 , 
- 1 , 



F 2 
1 

12 
1 



M P {s) 



1-4- 



Jp(s) - 



6 



+ 



288tt 2 



16tt 2 
1 

32^ 



, . , apis) — 1 
a P (s) In ^ : — +2 



<tp(s) + 1 



In ^ 

/i 2 



(B.4) 



Notice that the (3 parameter includes the contribution of the isospin breaking p — u mixing 
through Q pu) = -3.3 x l(T 3 GeV 2 [56], and F v and G v are defined in Eq. (21). Moreover 
the asymptotic constraint on the Nq — > oo vector form factor indicates FyGy ~ F 2 [54]. 
The mixing between the octet and singlet vector components employed in the construction 
of the 1 = component of the kaon vector form factors is defined by : 







cos 6y — sin ^ j I v 8 
sin^v cos6>y J \ v o 



(B.5) 



and we will use ideal mixing, i.e. 6y = 35°. 
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